Testing for a Super- Acceleration Phase of the Universe 
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We propose a method to probe the phenomenological nature of dark energy which makes no 
assumptions about the evolution of its energy density. We exemplify this method with a test for 
a super-acceleration phase of the universe i.e., a phase when the dark energy density grows as the 
universe expands. We show how such a phase can be detected by combining SNla (SNAP-like) 
and CMB (Planck) data without making any assumptions about the evolution of the dark energy 
equation of state, or about the value of the matter density parameter. 



I. INTRODUCTION 

Observations of distant Type la supernovae imply a 
presently accelerating universe 0,0- In the standard 
cosmological model this is accommodated by introducing 
"dark energy" , a component which interacts with the rest 
of the universe only gravitationally and has a significantly 
negative pressure which leads to the acceleration of the 
universe. 

In order to accelerate the universe, the dark energy 
component must have an energy density which decreases 
(if at all) much more slowly than matter density as the 
universe expands. Current data favor a dark energy den- 
sity which is almost constant or increasing with time 
[J IE El Jp and exciting results can be expected 
in the future |°l llOlllllll^| . We label the phase when the 
dark energy density is increasing with the expansion of 
the universe as super-acceleration. 

In this paper we discuss a method to probe phe- 
nomenological properties of dark energy without any as- 
sumptions about the evolution of dark energy density. 
We use this method to formulate a test to ascertain if 
the universe evolved through a super-acceleration phase 
by combining SNla and CMB experiments. We quote 
our results in terms of two variables w sn and w cm b which 
are the best fit constant dark energy equation of state 
parameters for the SNla and CMB experiments respec- 
tively. Our method does not rely on the equation of state 
parameter being constant. In fact, we do not need to as- 
sume anything about the evolution of the dark energy 
equation of state, or about the value of the matter den- 
sity parameter. 

We will outline our method in the context of the spe- 
cific example of testing for super-acceleration using SNla 
and CMB experiments. We will end with a discussion 
of how our method can be generalized to other questions 
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one might ask about dark energy using SNla, CMB and 
other experiments. Questions of the nature, "did the 
dark energy density increase with cosmic expansion?" 
are perhaps easier to answer than some questions pre- 
viously asked in the literature. And yet they have the 
potential to revolutionize our thinking about dark en- 
ergy. We stress that answering such questions does not 
necessarily require measuring the complete dark energy 
density evolution; the key point is that we are asking for 
less information. 



In a nutshell, to apply this method one needs to fit 
a constant equation of state w to each of SNla and 
CMB observations and plot conditional probabilities in 
the w versus J7m plane (e.g., Fig. 1). The probabili- 
ties are conditional in the sense that they are normalized 
at each value of £Im- If the upper contours cross be- 
low w = — 1 then super-acceleration has been detected. 
We stress that the power to make such robust state- 
ments comes from combining different measures (SNla 
and CMB) which have complementary dependencies on 
the dark energy equation of state and matter density. 



In Section [n] we start with a discussion of dark en- 
ergy and define super-acceleration rigorously. We then 
outline the general characteristics of models which could 
describe a super-acceleration phase. Given that dark en- 
ergy is a total mystery even if it never went through 
a super-acceleration phase, our point of view is that 
super-acceleration is a question to be observationally set- 
tled. It is satisfying to note that this will be possible 
with upcoming experiments. We explain our evolution- 
independent method in Section IIIII We then proceed 
to forecast the possibilities with regard to measuring a 
super-acceleration phase from future missions in Section 
FVl We end with a discussion of how the method can be 
generalized. 
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II. DARK ENERGY AND 
SUPER- ACCELERATION 

In this section we will consider dark energy in a gen- 
eral setting (without reference to a theory of gravity) 
and define it in terms of its phenomenological properties. 
This will allow us to motivate super-acceleration as a part 
of phenomenologically viable parameter space. Finally, 
we will consider some models where a super-acceleration 
phase might arise. These models are discussed with an 
aim to better understand the requirements that super- 
acceleration imposes. 



A. Generalized Equation-of-state 

The pressure of the dark energy (DE) px is parameter- 
ized in terms of the equation of state wx which relates 
it to the physical energy density px as px = wxpx- 
In terms of wx the evolution of the DE density with 
scale factor a is px{o) oc exp(— 3 J G?lnet'(l + wx(a')). 
If wx < —1, then the DE density is actually increasing 
as the universe expands, which seems counter-intuitive. 
In the context of a scalar field (4>) with canonical kinetic 
(0 2 /2 > 0) and potential (V(<f>) > 0) terms, w x < -1 
is unphysical. This is easily seen by rewriting wx as 
l + wx = 2<j) 2 /{<j) 2 + 2V(<j))} which is positive if V(4>) > 0. 

The SNIa observations provide us with a magnitude- 
redshift relationship. If the SNIa can be standardized as 
sources with known luminosity, then they measure the 
luminosity distance d^(z) in units of 1/Hq where Hq is 
the present expansion rate of the universe. Within the 
context of homogeneous and isotropic universes, we can 
then write 



s H d L (z) 



(1 + *) 



dz 
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(1) 



where h(z) — (a/a)/Ho is the expansion rate of the uni- 
verse normalized to unity today and z = 1/a — 1 is the 
redshift. Acceleration refers to the condition a > 0. We 
can parameterize the expansion rate as 



h 2 (z) = {1+zf 



9. 
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This is completely general and it does not tie us down 
to a theory of gravity (see [TJ EH EJ, EE El for ex- 
amples of studying cosmology in this spirit). Any h(z) 
can be written in terms of wx (z) as above. Eq. is a 
phenomenological fit; we need a theory (for gravity, dark 
matter and dark energy) to interpret 9m and wx{z). 

We can differentiate Eq. [21 to derive the acceleration 
of the scale factor. This yields 

2d/a = -(1 + 3w x (z))h 2 {z) + 3w x (z)(l + z) 3 9 M ■ (3) 

The universe is accelerating if w x (0) < — 1/[3(1 — 9m)}- 
Also, for -1/[3(1 - 9 M )] > wx(0) > -1/(1 - &m) wc 
have fi(O) > but h(0) < 0. 



B. Super- Acceleration 

We will refer to an expansion phase with w x (z) < — 1 
as super-acceleration j33|. Super-acceleration implies 
that the dark energy contribution to the expansion rate 
is increasing with time and this leads to interesting sce- 
narios for the ultimate fate of the universe [l7l Il8| . In 
particular, it is possible that the universe could end in 
a "Big- Rip", with the time scale for such catastrophic 
events set by I/Hq. In order to draw such conclusions, 
the future evolution of wx(z) must be known. 

We are interested in the issue of whether future obser- 
vations will be able to ascertain that the universe went 
through a super-acceleration phase. The main motiva- 
tion for thinking about super-acceleration is the simple 
fact about dark energy - we don't know what it is. Given 
our current understanding, even if we know what wx{z) 
exactly is, dark energy will still be a mystery. 

One of the troubling things about super-acceleration 
is that it seems to violate causality. The magnitude of 
the pressure is larger than the energy density and hence 
one might expect to transmit information faster than the 
speed of light. While this is a possibility, it is by no means 
guaranteed by having \px\ > Px = Sui^Hq [h 2 {z) — (1 + 
z) 3 9m}- wx{z) determines the global evolution of px', it 
does not tell us how small local perturbations propagate. 



C. Modeling Super- Acceleration 

Given the repercussions and the unintuitive nature of 
wx{z) < — 1 systems, it is of interest to investigate mod- 
els which give rise to a super-acceleration measurement. 

1. One possibility is that the universe is accelerating 
due to a cosmological constant (or perhaps some 
dynamical scalar field), but the effective wx at 
some epoch goes below -1 because of new physics 
we are unaware of. An example of new physics E3 
could be the coupling of photons to axions (not the 
QCD axion) [2(J which has been proposed to ex- 
plain the dimming of supernovae without the need 
for an accelerating universe. 

2. A second (troubling) possibility is that the effec- 
tive wx is less than -1 because of systematic ef- 
fects. This is similar to (1) except it is not new 
physics but some systematic effect that results in 
the "apparent" super-acceleration. 

3. wx < —1 is a possibility if one postulates that dark 
energy is a scalar field with non-canonical kinetic 
terms 0, El HE 113 ■ One must be careful about 
instabilities in the theory in this case [24| . 

4. Another class of models appeals to a modification 
of gravity on scales < 1/Hq to obtain w x (z) < — 1. 
Like the other models in this section, this class of 
models does not solve the "why now?" problem 
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since the scale Hq is put in by hand to make sure 
that gravity is only modified today. There are se- 
vere constraints on modified gravity theories in the 
early universe |14| . The main impetus for consider- 
ing modified gravity theories is provided by brane- 
world models wherein it might be possible to get 
super-acceleration [l^.l25| . 



III. DETECTING SUPER- ACCELERATION 

In this section we lay out the method by which fu- 
ture experiments can unambiguously detect a super- 
acceleration phase, independent of the functional form 
of w(z). First we show that if a constant equation of 
state is fitted to SNIa observations and a result less than 
— 1 is obtained then cither the true equation of state must 
at some point dip below —1, or the assumed matter den- 
sity is incorrect. We then show the same thing for CMB 
observations. Finally we show how SNIa and CMB in- 
formation together can overcome the need to assume a 
value for the matter density. Throughout we assume a 
flat universe. 



A. Information from Supernova observations 

Given a set of SNIa observations our method uses con- 
ventional x 2 fitting to obtain the best fitting constant 
equation of state for each Qm (see next section). The aim 
of the following paragraphs is to show the significance of 
this fitted constant equation of state, given that the ac- 
tual equation of state is of unknown functional form. 

A SNIa experiment measures the luminosity distance 
De out to a certain set of redshifts. As usual, we define: 



D E (z) = 



1 



dz 



Mm Jo + z) 3 + r E (z;n M ) 



(4) 



where r E (z;Q M ) — (^x/^m) exp[3 J* dz'(l + 
w x (z'))/(l + z')]. The subscript E in Eq. E| ex- 
plicitly denotes that this is the quantity measured 
by an experiment. We are further assuming that the 
experimental magnitude redshift relation is an unbiased 
estimator of the actual luminosity distance. 

The next step in our method is to fit the distance 
De{z) with a distance Dp(z, fijy, Aw) which assumes 
a constant equation-of-state, wx(z) — w sn = — 1 + Aw. 
Thus Dp is has the same form as De but with r F re- 
placed by r F (z; Q M ) = {Q x /^m)(1 + z) 3Aw . The best- 
fit Aw is found by maximizing the likelihood assuming 
gaussian errors (appropriate for SNAP) for the measured 
distances. 

We give a more rigorous derivation in the following 
paragraphs, but to gain a qualitative understanding first 
consider the weighting function approach of Saini et. al. 
[2(| . They showed that to a good approximation a con- 
stant fitted equation of state from supernovae observa- 
tions is simply a weighted average over redshift of the 



true equation of state, where the weighting function is 
always positive. From this it is clear that if the con- 
stant fitted equation of state is less than —1 then the 
true equation of state must also go below —1 for some 
redshift interval(s). The above is only true if the correct 
matter density has been assumed. 

We now have the task of computing the best-fit Aw for 
a given De{z) at a fixed Qm- This is found by minimizing 
the x 2 at fixed Qm with respect to Aw. The resulting 
equation does not have a simple solution. Our aim is to 
figure out when the best-fit Aw < 0. To achieve this, 
we replace Dp in the minimization equation by its a first 
order Taylor expansion in Aw, i.e., 

D F (z,n M ,Aw) ~ D A (z,Q M ) + Aw G A (z,n M ) . (5) 

D\ is defined the same way as Dp but with Aw = 0, 
i.e., D\(z,Hm) is the distance in ACDM cosmology. We 
define G F = dDp/dAw and the function G\(z,Qm) = 
G F (z, Q M , Aw = 0). For fl M > 0.1 and < z < 1.7, the 
error in the approximation given by Eq. [SJis smaller than 
3% for -0.5 < Aw < 0.5. We note that Saini et. al. HJ 
have shown that this kind of linear approximation works 
extremely well even when Aw is not constant (in which 
case we need to take a functional derivative with respect 
to Aw(z).) 

We note three facts. One, we are only approximating 
the fit (to ascertain whether Aw < 0) and not making 
any assumptions about the actual wx(z). 

Two, the error in approximating Dp according to Eq. 
[S] turns out to be unimportant. At the end of this sec- 
tion, we show explicitly that the only effect of this ap- 
proximation is to make our results (on detecting super- 
acceleration) somewhat conservative. 

Three, we could just as well use Dp = D\ + AwG\ 
as our exact fit and it would not change any of our argu- 
ments. We chose not to go this route because the above 
fit is unintuitive and has unphysical behavior for large 
Aw. 

Applying the approximation in Eq. [5] to the equation 
resulting from the minimization procedure, we get the 
following implicit equation for the best-fit Aw: 



Aw = 



Sj [D E (zj) - DaIz^Qm)} Gp{zj, n Ml Aw)/a 2 (zi) 
Gp(zi,Q M , Aw)GA(z i ,fl M )/o- 2 (z l ) 

(6) 

Both G F and Ga are negative for all z. What Eq. [fjjsays 
is that if the best fit Aw < 0, then D E {zi) > Da(zi, Qm) 
for a finite number of Zi or U,m is wrong. This in turn 
implies that if is the true value, then wx (z) must be 
less than -1 for some redshift range. 

A point of note here is that for each Qm we have a 
corresponding best-fit Aw(Qm)- We arc not minimizing 
with respect to both Qm and Aw. The f2jvf resulting 
from a joint minimization need not be the true matter 
density and in fact it could be that the best-fit Aw < — 1 
even though the underlying wx(z) never enters a super- 
acceleration phase [2^ • 
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FIG. 1: The shading shows the two-dimensional probability 
distribution from SNAP-type supernova data in the w an -flM 
plane. The solid lines show 68 per cent limits at each Qm 
from the same simulated data. The dashed lines show 68 
per cent limits at each Qm from a Planck-type simulation, 
marginalized over 5 other cosmological parameters. The input 
model has Qm = 0.3 and w(z) — —1 (a constant). 



We expect a few hundred supernovae from SNFactory 
|2m and the low redshift Carnegie Supernova Program 
|29| and that will enable calibration of the luminosity 
distances to around 0.5 per cent. If the ~6000 supernovae 
that should be seen by GAIA |30| were all followed up 
then this could reduce the uncertainty to 0.2 per cent. 

This additional uncertainty will reduce our ability to 
determine if super-acceleration occurred. We take the 
conservative (and not fully optimized) approach of as- 
suming the worst possible value for the calibration pa- 
rameter and feeding it through the analysis described 
above. If the luminosity distance calibration were actu- 
ally a factor of 1.015 (3x0.5 per cent) higher than esti- 
mated from the low redshift sample, then the supernova 
contours would shift from the solid lines in Fig. [2] to the 
dotted lines. This factor corresponds to the true calibra- 
tion factor being about three <r away from that estimated 
by the upcoming experiments. If there were significantly 
more supernovae followed up, then this would improve 
the picture. For example, we get the dashed contour 
with a factor of 1.006. 



How are our results affected by the small error in ap- 
proximating the fit luminosity distance fljr by a first or- 
der Taylor expansion? To answer this question let us 
write 

D F (z,Q M ,Aw) = {D A (z,n M ) + AwG A (z,n M )) 

x (l-e(z,Q M ,Aw)). (7) 

For the interesting parameter space, e(z, Qm, Aw) is 
smaller than about 0.03. Further, for all practical pur- 
poses, this approximation overestimates Dp. One can 
then show that the effect of this is to make Aw smaller. 
Thus the conclusions drawn from Eq. are conserva- 
tive in the sense that Aw is more negative than what we 
would find from Eq. 

An interesting exercise is to fit the supernova distances 
with Dp(z, Qm, Aw) = D\(z,Q M ) + Aw G\{z, Qm) ex- 
actly. We can then run through the same arguments as 
before and arrive at the same conclusion. In particular, 
we computed the new contours with w sn = — 1 + Aw and 
the quantitative results changed very little as expected. 

B. Low redshift supernova calibration 

In the discussion above, we have assumed that we can 
measure luminosity distance Dj,asa function of redshift 
from supernovae. However in practice we can only mea- 
sure magnitudes m 

m = M + 5logD L (z,Q m ,w(z)), (8) 

where the calibration factor M. is unknown unless we 
have a very large number of low redshift supernovae. 



C. Information from CMB observations 

The CMB is sensitive to the dark energy mainly 
through its effect on the angular diameter distance to 
the CMB, and therefore the positions of the acoustic 
peaks 31] . It also has an effect on the largest scales 
since dark energy changes the Integrated Sachs Wolfe ef- 
fect, through both the change in the expansion rate and 
clustering. Here we propose to ignore the CMB infor- 
mation on large scales so that there is no need to make 
assumptions about the clustering of dark energy. 

In the case of the CMB, the approximation given in 
Eq. [3] does not work well. However we do not need 
the approximation since the CMB measures the angular 
diameter distance out to just one redshift, that of the 
last scattering surface. The equation to solve for the 
CMB is Dp(oc, Qm, Aw) — De(oo) (neglecting the small 
error on Dp)- Although the relation between Aw and 
De{z) is non-trivial here, the argument is simple. For the 
true Qm if the relation Dp(z, Qm, Aw) — Dp(z) implies 
that Wcmb = — 1 + Aw < — 1 then a super-acceleration 
phase must happen over a finite range of redshift. This 
is simply a consequence of the fact that Dp(z, Qm, Aw) 
is a monotonically decreasing function of Aw. 

Qualitatively, when a constant equation of state is fit- 
ted to CMB data it measures roughly a weighted average 
of the true wx (z) where the weighting function is always 
positive [2^|. Therefore if the fitted w cm b gives a value 
less than — 1 then the true equation of state must itself 
dip below —1. As for the SNla argument, this is only 
true if the value of Qm has been guessed correctly. 
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FIG. 2: The contours in this figure are computed the same 
way as the contours in Fig. The solid supernova contour 
assumes that the calibration error is negligible. The dashed 
lines are 68 % limits assuming a worst case calibration pa- 
rameter estimated from 6000 low redshift supernovae. The 
dotted line is the same as the dashed but for 300 low redshift 
supernovae. 



D. Combination 

As discussed above, using SNIa alone it is not possible 
to uncover the nature of dark energy simultaneously with 
estimating the matter content Qm- This is also true for 
CMB data alone. Therefore we need to combine informa- 
tion from SNIa and CMB. Given observed data we can 
determine the limits on w at each Qm for both SNIa and 
CMB and plot this conditional likelihood in the w-Qm 
plane. In general the best-fit w at some Qm derived from 
SNIa will be different from that derived from CMB since 
they probe the equation of state at different redshifts. 
For illustration we use contours which encapsulate 68% 
of the likelihood throughout, however the method can be 
carried out for whatever confidence level is required. 

The key point of our argument is that if, for every value 
of Qm, the contour for either u) cm b or w sn falls below -1, 
then it must be that dark energy went through a phase of 
super-acceleration for a finite period of time. We explain 
how this follows from the previous sections by discussing 
two examples. 

Consider the example in Fig. ^ Suppose that we do 
not know the true value of Qm- We get around this 
problem by considering each possible value for Qm in 
turn. For each Qm value we assess whether, if this were 
indeed the true value of Um, we would conclude that the 
universe underwent super-acceleration. If for all values 
of Qm we reach the same conclusion, then we can give a 
definitive answer, independent of the value of Qm- 

As an example let us first consider the value Qm = 0.2. 
In this case w sn > — 1 and w cm h < — 1. If this were the 
true value of Q m then we would have to conclude that at 



some point in the history of the universe, w(z) < —1. (In 
fact since the CMB probes the value of w(z) at higher 
redshifts than SNIa then this would imply that w was 
once less than —1 and has since risen above —1, but these 
details arc not important for our argument.) We would 
reach this same conclusion for all matter densities with 
Qm ^ 0.28. Now consider the value Qm — 0.4. If this 
were the true value of Qm then we would again conclude 
that for some redshifts w(z) < —1, since u) sn < — 1. We 
would reach the same conclusion for Qm ~ 0.305. How- 
ever, for Qm — 0.3 then both data sets allow w > — 1, 
therefore we cannot be sure whether or not there was 
super-acceleration. Our conclusion from Fig.^is that ei- 
ther 0.28 ;$ Qm ~ 0.305 or there was super-acceleration. 

Now consider future observations providing the solid 
(SNIa) and dashed (CMB) lines in Fig. g] Running 
through the same arguments as above, we conclude that, 
no matter what the value of Qm, there must have been 
super-acceleration. This can be seen at a glance by the 
fact that the upper contours cross below —1. 

The method we have just outlined determines if the 
universe underwent a super-acceleration phase without 
any assumptions about the evolution of dark energy den- 
sity. It is clear from the above discussion that one cannot 
test for a super-acceleration phase if the deviation from 
wx = — 1 is arbitrarily tiny. How well can we do with 
future experiments? We turn to this question next. 



IV. FORECASTS 

In order to make a prediction about how well future 
experiments will do, we need to generate mock data (De) 
which in turn requires that we make some assumptions 
about the evolution of wx(z). This is not a major hur- 
dle (as we show later), mostly because we interested in 
getting a rough sense of what future experiments can 
achieve. In the next section, we will discuss the proce- 
dure used to generate mock data assuming a particular 
wx{z) function. We generate data relevant for SNAP- 
like and Planck experiments. 



A. SNIa data 

We assume that the luminosity distance will be mea- 
sured to an accuracy of 1.4 per cent in each of 50 red- 
shift bins evenly spaced between redshift 0.1 and 2. First 
we simulate data from a model with w(z) — —1 and 
Q M = 0.3. 

On fitting the two parameters, w sn and Qm we obtain 
a tight constraint on the two parameters indicated by the 
shading in Fig. ^ 

For our method the quantity of interest is not the two- 
dimensional probability (as shown by the shaded con- 
tour) but the confidence limits on w sn calculated for a 
given Qm- These limits are shown by the solid lines on 
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Fig. ^ an d as expected the contours extend to w sn < — 1 
in the direction of increasing Qm- 



B. CMB data 

We generate CMB temperature anisotropy (TT) data 
using the parameters for the Planck HFI 353 GHz chan- 
nel: sensitivity AT /T = 14.4 x 10~ 6 , solid angle per pixel 
47r/5.94 x 10 6 (where 5.94 x 10 6 is the number of pixels 
on the sky) and 5 arcmin FWHM beam and assume 100 
per cent of the sky can be used. To speed up the com- 
putation we bin the observations with A£ = 10 and go 
up to £ max — 2500. We also marginalize over the Hub- 
ble constant, the baryon content, and the amplitude and 
spectral index of the primordial scalar fluctuations. We 
ignore tensor modes and do not use CMB polarization 
information, although in principle these could both be 
included in the analysis. 

The CMB data essentially constrain the angular diam- 
eter distance to last scattering z ~ 1100 and therefore 
there is a perfect degeneracy between u) cm b and Q~m- In- 
cluding the large scale ISW effect (the decay of gravita- 
tional potential due to dark energy domination) or lens- 
ing on small scales would break this degeneracy to some 
extent. However, a reliable estimation of the degeneracy 
breaking would require some assumptions about the clus- 
tering of dark energy. In light of the uncertain nature of 
the models of dark energy with wx{z) < —1, we have 
decided to not use the lensing signal or the ^ < 50 power 
spectrum (where the ISW effect is important). We again 
plot contours at the 68 per cent confidence limits for each 
CIm value, although the more conventional contours at 
constant two-dimensional probability look similar. 



C. Requirements on w(z) for detection 

In the simplest approximation the shape of the con- 
tours in Fig.QJis independent of the shape of wx(z). For 
now we assume this to be true and assume the shape is 
the same as that for a constant wx- We discuss the va- 
lidity of this approximation in the following subsection. 

We produce plots similar to Fig. ^ but simulate data 
for models with various values for the constant equation 
of state. We fix fl^ = 0.3 throughout and allow differ- 
ent values for the equation of state for the CMB simu- 
lations and the SNIa simulations. This corresponds to 
the case where the equation of state is varying with red- 
shift and therefore the effective value is different for the 
two probes. For each assumed u> cm b we calculate the 
w sn for which the upper CMB and SNIa 68% contours 
cross below —1, and therefore it is possible to detect 
super-acceleration. This is shown by the solid line in [21 
For this value of w sn and lower it is possible to detect 
super-acceleration. The dashed and dotted lines show 
the equivalent for Q~m = 0.2 and Q~m = 0.4. 




FIG. 3: The value of w sn below which it is possible to detect 
super-acceleration, for a given w cm b and Qm- The solid curve 
corresponds to a fiducial matter density of Qm = 0.3, the 
dashed line corresponds to Qm = 0.2 and the dotted line 
corresponds to Qm = 0.4. The cross shows the w a n at which 
super-acceleration can be detected for an evolving equation 
of state with a jump in w x (z) at redshift 0.3 for Qm — 0.3. 



This takes into account the error bars predicted for 
SNAP-like and Planck-like experiments. For example, if 
the true value of Q~m is 0.3 then if we measure w} C mb = — 1 
then w sn must be less than —1.3 if we are to detect super- 
acceleration from these future experiments. Note that it 
is not necessary for us to assume a value for VIm in order 
to detect super-acceleration; however the true value of 
Q,m will affect how easy it would be to detect such a 
phase. 

In making Fig. 01 we have assumed that the calibration 
error is negligible. Depending on the calibration error, 
the curves in Fig. [21 will shift down somewhat. The mag- 
nitude of this effect may be ascertained by looking at the 
shift of the contours in Fig. [21 



D. Insensitivity to assumed fiducial wx 

To calculate the lines in Fig. [21 we assumed that the 
shape of the contours (e.g., of Fig. ^) are independent 
of the form of w(z). We take two approaches to assess 
the reliability of this approximation. First we calculate 
the mathematical requirements for this approximation to 
be true. Then we take an extreme example of a sharply 
changing w{z) and show how even in this case, our pre- 
dictions in Fig. [21 are not altered much. 

To determine the mathematical requirements for our 
approximation let us assume we simulate mock data with 
a given value of matter density denoted by VL M ("true 
matter density" ) and a particular evolution for the dark 
energy density. For this De{z) we find the best fit 
Aw(Um) (as described above in Section ITTTf) . We are 
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0.45 



FIG. 4: (a) The equation of state as a function of redshift for an extreme strongly evolving model. A large jump occurs at 
2 = 0.3. (b) Solid blue lines: 68 per cent confidence limits on w an for each Qm value from a SNAP-type simulation with w(z) 
given by Fig. 3(a). Dash-dot lines: the simulation has w(z) — —2.38, the best fitting constant equation of state model to the 
supernova data. The solid and dash-dot lines are quite similar, despite the fact that the wx evolution is so different in the 
two cases. Since this w(z) was devised to give the same angular diameter distance to the CMB as a w — — 1 model, the CMB 
contours (dashed lines) are unchanged from Fig. 1. 



interested in the question of how Aw(Qm) changes due 
to a change in De{z) with the constraint that Aw(0' M ) 
is unchanged. Let us consider the following change 
D E (x) -> D E {x) + SD E (x). Then 



Aw(Q M ) = Aw(Q M ) + 

Si SD E {zi)G F (zi, nyWgfe, Qm, Q M )I^{^) 
£j G F (zi, Q M )G A (zi, 9, M )/cr 2 {zi) 



,(9) 



where we have dropped the explicit dependence on Aw 
in the right hand side of equation. In the above Eq. 03 



w F (z,n M 

,Qm 



,Aw) 



G F {z, Qm, Aw) 
G F {z, Q^, Aw) 



SI* 



1) 



din G F (z, Qm, Aw) . 



d\nQ 



M 



(10) 



where the last approximation is valid for small depar- 
tures from Q l M - This is in fact the situation we will be 
interested in. In this regime, W f (z,Qm,Qm) is sensibly 
constant. If W F (z) is truly constant, then Aw(Qm) = 
Aw(Qm)- Thus we have shown that a Aw — Qm con- 
tour derived from a particular D E (z) is not significantly 
changed due to changes in D E (x) which leave the best- fit 
Aw(Q M ) unchanged. 

We can illustrate the above formal derivation with an 
example. Consider the w x (z) given by the solid line in 
Fig. ||| (a), which has w x (z) = -2 for < z < 0.5, 
Wx{z) = for 0.5 < z < 9 and wx{z) = —1 at higher 
redshifts. We keep Q m — 0.3. This equation of state 
gives the same angular diameter distance to z — 1088 
as the wx(z) — —1 simulation, so the CMB constraints 
are unchanged. Since the supernovae are most sensitive 



to low redshifts, the effective equation of state seen by 
the supernovae is closer to -2. The resulting conditional 
w versus Qm plot is shown in Fig. Ub) . The solid lines 
show the constraints from SNAP-like observations from 
a simulation using the w(z) shown in Fig. Ufa). If we 
instead simulate SNAP-type data from a model with a 
constant equation of state, fixed to the best fitting w sn 
at Qm = 0.3 we obtain the dot-dashed line. By design 
the two coincide at Qm = 0.3, and overall they are very 
similar, as expected if our approximation is good. 

We now modify this jumpy wx{z) model to find the 
limiting case where the upper contours only just cross 
at w = —1. We keep the transition redshifts the same, 
and vary the high and low values of w x , maintaining 
the condition that u) C mb = — 1 (at Qm = 0.3). We find 
that a model with wx(z) = —1.1 (0 < z < 0.3), and 
wx{z) — —0.9 (0.3 < z < 9) meets this criterion. This 
has w sn = —1.03 and and so for this shape of wx{z) 
the solid line of Fig. |3| should actually pass through the 
small cross. This illustrates the effect of different shapes 
of w(z) on the result in Fig. O 



V. CONCLUSIONS 

Throughout we have concentrated on the question of 
super-acceleration. However our method can be extended 
to other similar questions. For example it can answer the 
question of whether or not w(z) ever went above —1. If 
the lower contours on the conditional w versus Qm plot 
cross above —1 then we can be sure that at some point 
w(z) was greater than —1. Alternatively we can answer 
questions about different values of w. Denoting the w 
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value where the contours cross as w u (upper contours) 
and wi (lower contours) we can say that at some point in 
the history of the universe w(z) < w u , and at some point 
in the history of the universe (not necessarily the same 
point!) w(z) > w\. In fact the line in the w, £Im plane 
need not even be horizontal. For example, one could 
imagine testing for whether w < — 1/[3(1 — £1m)]- 

It is possible that this method could be extended to 
other methods of constraining the equation of state of 
dark energy, for example cosmic shear and cluster counts. 
This might make it easier to detect a super-acceleration 
phase. 

In conclusion, we have discussed a method which can 
detect the super-acceleration of the universe without any 
assumptions about the evolution of the dark energy den- 
sity, or about the present matter density. This exam- 
ple makes it clear that despite the intrinsic degeneracies 
in supernova and CMB observations and our ignorance 



about the evolution of dark energy density, it will be pos- 
sible to make fundamental and robust discoveries about 
the phenomenological nature of dark energy with upcom- 
ing SNAP-like and Planck observations. 
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